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Abstract 

Current problems with the solar model may be alleviated if a significant amount of dark matter 
from the galactic halo is captured in the Sun. We discuss the capture process in the case where the 
dark matter is a Dirac fermion and the background halo consists of equal amounts of dark matter 
and anti-dark matter. By considering the case where dark matter and anti-dark matter have 
different cross sections on solar nuclei as well as the case where the capture process is considered to 
be a Poisson process, we find that a significant asymmetry between the captured dark particles and 
anti-particles is possible even for an annihilation cross section in the range expected for thermal 
relic dark matter. Since the captured number of particles are competitive with asymmetric dark 
matter models in a large range of parameter space, one may expect solar physics to be altered 
by the capture of Dirac dark matter. It is thus possible that solutions to the solar composition 
problem may be searched for in these type of models. 
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I. INTRODUCTION 


Over the past few decades, a large amount of evidence in support of the existence of dark 
matter (DM) has been assembled, see e.g., refs IH |2]. To describe the properties of DM 
has become one of the main issues not only in cosmology but also in particle physics, since 
the Standard Model leaves no room for its existence, which means extending it is necessary. 
Some of the most studied particle candidates of DM are weakly interacting massive particles 
such as the lightest neutral supersymmetric particle, the neutralino (an in-depth review can 
be found in ref [3]). 

If DM interacts with regular matter, it may scatter in astrophysical bodies such as the 
Sun and become gravitationally bound. Over time, this can lead to a large accumulation of 
DM, the effect of which used to be one of the proposed explanations of the solar neutrino 
problem [1]. This idea was later discarded in favour of neutrino flavour conversion by the 
neutral current phase of the SNO experiment [5]. Since the downward revision of heavy 
elements in the standard solar model |6j, theoretical predictions and observations of helio¬ 
seismology do not match and the Sun now faces the solar composition problem [7j. Again, a 
proposed solution is DM trapped in the Sun due to its interaction with regular matter. The 
idea is that DM particles collide with nuclei in the Sun, lose enough energy to become grav¬ 
itationally bound and to eventually settle in the solar core after additional scattering. Once 
in the core, DM can scatter off of the thermal distribution of nuclei and gain some energy 
which is then lost by scattering in the outer regions of the Sun, the result of which is a lower 
core temperature [8]. This would ultimately affect helioseismology and possibly provide a 
cure for the solar composition problem. The shift in temperature of the solar core would also 
change the solar neutrino fluxes that are observed by experiments, in particular the flux 
which varies as [S] . The effects of DM captured by the Sun has been extensively studied 
for various DM models [T0HT6] and also in other stars [IMI]. Annihilating DM can also 
provide a channel for indirect detection of DM as a flux of high energy neutrinos from these 
annihilations might be detectable on Earth. Such signals are being searched for by neutrino 
telescopes [22l [23]. In the case where the DM has self-interactions, DM particles already 
captured by the Sun provide a possibility of DM self-capture, which could lead to higher 
concentrations of DM inside the Sun. If DM annihilations produce high-energy neutrinos. 
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this would increase the number of annihilations and thus the signal in neutrino telescopes 
could be enhanced [21]. Energy injection into the Sun by the DM annihilation products has 
been studied, but requires the halo density to be many orders of magnitude higher than the 
local density for the star to be affected [25]. 


The models which achieve the highest numbers of captured DM particles in the Sun are 
models in which annihilation is suppressed or does not occur at all, the latter of which is the 
case in asymmetric DM (ADM) models, reviews of which can be found in [26l |2Z] . It was 
proposed in [T2] that self-interacting ADM may alleviate the solar composition problem but 
it was shown in osum that ADM models are incapable of doing so with the cross-sections 
allowed by direct detection experiments. Recently, there have been improvements of the 
standard solar model using DM has been achieved with long range interactions [2H] and 
momentum-dependent cross sections 123. Common for the models described above is that 
they contain only one species of DM particles. In this paper we will study the accumulation 
of DM in the Sun under the assumption that it is a self-interacting Dirac fermion and that 
the DM halo is composed of equal amounts of DM and anti-DM particles. This opens up a 
possibility for an asymmetry in the captured number of DM and anti-DM to occur, either 
due to different capture cross sections for DM and anti-DM on solar nuclei or due to random 
fluctuations in the capture process. Moreover, the more abundant species of captured DM 
can not fully annihilate with its counterpart which allows for large number of captured 
particles even for an annihilation cross section expected for thermal relic DM. If the size of 
the population of captured DM of the models considered here can compete with the ADM 
models, DM of the standard freeze-out scenario may alter solar observables. 


The remainder of this paper is organised as follows: In section |TT| we present the frame¬ 
work that governs the amount of DM and anti-DM present in the Sun and discuss the 
behaviour of the DM to anti-DM asymmetry in the limiting cases. Then, in section HI we 


apply our framework to the case of the Sun and in section IV we present our results. Finally, 
in section |V| we summarise our results and present our conclusions. 
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II. ACCUMULATION AND THE ASYMMETRY 


We model the total number of DM particles in the Sun N and the total number of 
anti-DM particles in the Sun N using the coupled system of hrst order differential equations 

N = c + CN + CN-TNN, (1) 

N = c + CN + CN - TNN. (2) 

Here, c and c are the capture rates on solar nuclei for DM and for anti-DM, respectively, 
while C and C will depend on the self-capture rates of DM and anti-DM and on the ejection 
rates of DM and anti-DM that is already trapped by incoming particles from the DM halo. 
The DM-anti-DM annihilation rate is given by T. In principle, C and C also include 
evaporation effects, but it has been shown that this is negligible for DM particles with a 
mass > 5 GeV pUl IHT] and we will neglect it in the following discussion. 

We dehne the asymmetry A as the difference between the DM and anti-DM numbers, 
A = N — N. It evolves according to 

A = d + DA, (3) 

where we have dehned d = c — c and D = C — C. Assuming negligible amounts of both DM 
and anti-DM at the birth of the Sun, A(0) = 0, we hnd 

A = f (4) 

Jo 

An important note is that the absolute value of A also represents the minimum amount of 
dark particles (DM or anti-DM) present in the Sun. The sum of DM and anti-DM, N + N, 
will at any time always be larger or equal to the asymmetry. 

Another important note is that there is a geometric limit for the self-capture rates. When 
N and N are large enough, our approach is no longer valid. When DM and anti-DM is 
trapped inside the Sun, they will accumulate in the center roughly inside a sphere of radius 
(derived in the next section). When + a^^N is larger than the cross-sectional area of 
the sphere, vrr^, every incoming DM particle will scatter off either trapped DM or anti-DM. 
When the distributions of trapped DM and anti-DM are identical, the fraction of collisions 
of DM on DM is = cr^^N/{cr^^N + while the fraction of collisions on anti-DM 

is given by + a^^N). The same argument holds for incoming anti-DM 
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except the total cross section can not exceed vrr^. In this case, the fraction of 

collisions on DM is given by with all N and N interchanged. Similarly, the fraction of 
collisions on anti-DM is given by again with N and N interchanged. When this occurs, 
the equation for A will fail since eqs. ([^ and (|^ must be corrected in order to take the 
geometric limit into account. 


II. 1. Intrinsic different capture rates 


When the scattering cross section of DM and anti-DM on solar nuclei are different, c and 
c are different and thus d ^ 0. We choose c to be larger than c (identifying anti-DM as the 
species with lower capture probability on solar matter). The solution for A is 

A = ;^(e“ - 1). (6) 


which has three interesting limits: 


dt if \Dt\ ^ 1 


A 


< if \Dt\ > 1,D < 0 


if|Dt|>l,D>0 


( 6 ) 


When Dt is small, self-capture is negligible and the asymmetry will be proportional to 
the difference in the capture rates. When D < 0, the system eventually stabilizes, since 
the additional capture of anti-DM by already captured DM balances the difference in the 
capture rate on normal matter. For D > 0, DM captures itself at a larger rate than anti-DM. 
Once this process becomes dominant, it leads to an exponential increase in the amount of 
DM captured in the Sun. 


II.2. Stochastically induced difference 

When the capture rates are equal (c = c), the amount of DM might at some point be 
larger than the amount of anti-DM simply due to random variations in the capture process 
which can be modelled by adding a white noise signal 6c to the capture rates 

c = co + 6c{t), {6c{t)) = 0, {6c{t)6c{T)) = s6{t -t). (7) 
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The white noise is normalized such that the expected number of captured DM particles and 
its variation matches those of a Poisson distribution, i.e., (n^) — (n)^ = (n). We hnd that 


(n) 

(n^) 


Co + 5c{r)dT ) = Cot 


Co + 5c{T)dT ) = Co t + st 


( 8 ) 

(9) 


and hence s = Cq. Using the same argument for the capture rate of anti-DM with a white 
noise signal 6c, we hnd 

d = 6d{t) = 6c{t) - 6c{t). (10) 


Since 6c and 6c are independent, 6d has the properties 


iSd{t)) = 0, iSd{t)S4(T)) = 2caS(t - t) 


( 11 ) 


The expectation value of the asymmetry A is zero, which should be expected since the 
probability of having an over-abundance of DM to anti-DM must be the same as that of 
having an over-abundance of anti-DM due to symmetry. To estimate the typical magnitude 
of the asymmetry, we can study the standard deviation of the stochastic variable A, given 
by A = (A^). We hnd that 

A^ = / [ - 1). (12) 

Jo Jo 

Thus, the limiting behaviour is similar to the case for intrinsic diherent capture rates 

if \Dt\ < 1 

if |Dt| > 1,D < 0 • (13) 

if \Dt\ 3> 1, D > 0 

The major diherence is that the short time limit \Dt\ 1 gives an asymmetry that is 
expected to grow with the square root of t rather than linearly and that the coefficients are 
related to Cq and D by a square root. For small Dt, we note that this result is precisely 
what would be expected from the diherence between two Poisson distributions of expectation 
value Cot while an equilibrium or an exponential growth occur for strong self-capture. 


A = 


y/'Jcot 

/_cg 
V D 


/Wf,Dt 
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III. DM AND ANTI-DM SELF-CAPTURE RATES 


For a self-interacting model of DM, the total capture of DM in the Sun will be the sum 
of a capture rate due to interactions with solar nuclei, a term proportional to the already 
captured DM and a similar term proportional to the number of captured anti-DM particles. 
The capture of anti-DM is completely analogous to DM capture although the rates may 
differ depending on the various scattering cross sections. The specihc formulas to compute 
the capture rates and ejection rates are presented in the appendix but the complexity of 
the capture rates C and C requires some discussion. In what follows, we assume that the 
time it takes for a captured DM particle to fall into thermal equilibrium in the solar core is 
negligible. 

Generally, the formula for the capture rates of DM by DM and anti-DM as well as the 
ejection rates are given by: 


1‘Rp) 


C{a) = / 47rr^ 


/(«) 


wVL du dr. 


u 


(14) 


The factor D is the rate at which a particle with velocity w will scatter at radius r and 
contain information on the probability that the incoming particle is captured and whether 
the target particle is ejected or not. 

The rate of capture of an incoming particle without ejecting the target particle is given 


by 

Cs{(t)= / 47rr^ / - -an{r){vlg^ —u^) du dr. (15) 

•JO Jo 

Depending on the collision, the scattering cross section a is either or The velocity 
of the particle before it falls into the gravitational potential is u and the velocity at radius 
r is then w = + Ugsc where Vesc is the escape velocity. 

The Knudsen number is a measure of the distance DM particles travel on average between 
collisions and is given by 

(Tini{r)^ . (16) 

The parameter is a length scale that describes the size of the distribution. For iF 3> 1, 
particles travel a great distance between collisions and an isothermal assumption of the 
distribution is justihed: 

^^iso(’^) = J^iso(0)e"'^^'’^/'''^, (17) 
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where niso(O) is the normalization constant, T is the temperature of the distribution and 
(j){r) is the gravitational potential energy at radius r from the core which is calculated from 



' Gm,M(r') 


( 18 ) 


Jo 


Here, M(r) is the total mass inside the sphere of radius r around the solar core; 


pr 

M{r) = / 47rp(r') 


dr'. 


(19) 


In the isothermal case, the normalization constant is uiso(O) = vr 2 r^^N and the length scale 


is defined as = ?>kTc/2'KGpcm^ assuming a constant temperature T(r) = Tc and density 
p(r) = Pc with Pc and Tc being the density and temperature in the solar center. Indeed, for 
a DM particle with mass = 5 GeV, the length scale ~ 0.05 Rq which shows that the 
vast majority of captured DM particles are concentrated in a very small volume in the center 
of the Sun. In the case of large scattering cross sections, R' -C 1, the particles scatter so often 
that they will be in local thermal equilibrium with the surrounding nuclei. The distribution 
for this case was derived in |8] which introduces a dependence on the temperature gradient. 
However, for the cross sections considered in this paper {asn < 10“^^ cm^), R' > 90 and 
so we will work with the isothermal distribution. Dehning e(r) = n{r)/N, the distribution 
can be written n{r) = e{r)N. For anti-DM, the radial distribution is taken to be the same 
except N N. 

The ejection rate of DM captured by in the Sun by collisions with DM or anti-DM from 
the halo is given by 



( 20 ) 


We must also take into account that while a dark matter particle is being ejected from the 
Sun, it is also likely that the particle from the halo is captured by the same process. The rate 
for this exchange occurring is Gexch = Geject — Geject 2 , where Geject 2 is the rate of ejections in 
which both the incoming halo particle and the particle from the Sun are ejected, given by 



( 21 ) 


The possible cases depending on velocity for capture and ejection are shown schematically 


in Fig. [T] along with the velocity distribution of the standard halo model. 
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FIG. 1. The behaviour of the functions wil. for different capture scenarios. The horizontal line 
represents the escape velocity Vesc in the center of the Sun, where most of the DM resides. The 
vertical size of the colored regions represents the quantity which needs to be integrated along with 
the distribution f{u)/u in order to yield the capture rates. The black curve shows shape of f{u)/u 
at Vq = 220 km/s for reference. 

As can be seen from this figure, the escape velocity will generally be so large that we 
expect that the self-capture will be dominant with some contribution from halo particles 
being captured while ejecting the target particle. 

Summarizing this, the capture rates C and C relevant for the evolution of the dark matter 
and anti-dark matter numbers in the Sun can be written as 


C C'eject(c’'xx) C*eject 2(f^xx)’ 

C = C's(^xx) A C*exch(c’'xx)■ 


( 22 ) 

(23) 


Here, the single self-ejection events do not appear in G, as the net change in the DM number 
in the Sun is zero for these events, but the full ejection induced by the opposite species from 
the halo must be taken into account as the capture of the halo particle does not compensate 
for the ejected one. For the capture of DM on anti-DM, the relevant quantities are the 
capture without ejection and the ejection of the target particle while capturing the halo 
particle. 
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The annihilation rate is compnted as [32] 


cRq 


r = (av) / 47 rr^e(r)^ dr, 


(24) 


where (av) is the thermally averaged annihilation cross section. As long as i? 03 >r^, the 
npper limit of the integral can be set to oo rather than Rq and the annihilation rate evalnates 
to 


r = (av) 


(27r)^Rr^ 


(25) 


IV. RESULTS 


In the following, we will make some explicit assnmptions in order to estimate the effects 
described in the previons sections. The velocity distribntion f{u) of the halo is assnmed to 
be a standard Maxwell-Boltzmann distribntion shifted to the solar frame moving throngh 
the halo at Vq = 220 km/s. It can be expressed as [1] 


/x(«) 


= 


u 




(26) 


It is assnmed that the DM and anti-DM components in the halo are identical and that the 
density of each are eqnal. They will then each have a density Px = Px ~ GeV cm“^, 
which is eqnal to half the total local DM density of 0.3 GeV cm“^ [531154] . The nnmber 
density of DM and anti-DM is therefore ny_ = p^/m^. 

The dark matter is assnmed to scatter with regnlar matter with velocity independent spin- 
independent (SI) and/or spin-dependent (SD) cross sections throngh effective operators. 
For the case of SD captnre in the Snn, we are interested in the bonnds on the SD DM- 
proton cross section. Limits on these cross sections have been set in varions direct detection 
experiments [331140] . In the DM mass range 10 — 1000 GeV, the limit on the SI cross section 
is asi ^ 10“^^ cm^ |35|. For smaller DM masses, the limits on the SI cross section weakens 
significantly. For a 5 GeV DM particle, asi ^ 10“^° cm^. The limits on the SD cross section 
in the mass range 10 — 1000 GeV is asD ^ 10“^® cm^ [37|. For a 5 GeV particle, the bonnd 
is slightly rednced to asn ^ 10 “^^ cm^ [40] . 

Limits on the self-interaction of DM comes from astrophysical sonrces. When galaxy 
clnsters collide, drag forces acting on the gas while the DM passes throngh nnhindered 
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FIG. 2. The self-capture rate Cgjo (black line), the exchange rate Cexch/o" (red line) and the 
ejection rate Geject 2/<7 (blue line) as a fnnction of mass. 

would produce an offset in the mass and gas distribution of the clusters, the size of which 
can be used to put upper limits on the self-interaction of DM. In [H], one such collision 
was analysed and and set an upper limit on the self-interacting cross section of o^xxl^x ~ 
2 ■ 10-24 cm2/GeV. 

The relic abundance of DM has been precisely derived from WMAP [12] and Planck ^3] 
experimental data. The thermally averaged annihilation cross section can be related to this 
relic abundance by solving the Boltzmann equation which is done in e.g., ref |3] and is here 
taken to be 3 ■ cm^/s. 

As a model of the Sun, the AGSS09 solar model m is chosen. It contains the mass and 
radial distribution of elements up to Ni. The solar age is taken to be t© = 4.5 byrs. 

The capture rate of a DM particle with mass rrix = 5 GeV is calculated to be at most 
lO^’^ s”^ for a SD cross section of 10“^^ cm^ and 2.9 ■ s”^ for a SI cross section of 

10-4® cm2. Pqj, ^ p)]y|- of mass 10 GeV, the bounds push SI capture down by four 

orders of magnitude even though for a hxed SI cross section, the capture rate is only reduced 
at the percent level. Thus, the SD cross section allows for higher capture rates for all masses 
in the range 5-1000 GeV. 

Fig. shows the values of Cg, Gexch and Geject 2 as a function of DM masses between 
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5 and 1000 GeV. While Gexch is roughly 20 times smaller than Cg, C'eject 2 is almost 15 
orders of magnitude lower. This is not surprising as the escape velocity is very large where 
the DM resides and particles with a large velocity in the halo are exponentially suppressed 
(cf. hg.g. In the case of the Sun, Geject 2 may therefore be neglected, since Cg and/or Gexch 
will be completely dominant depending on the relative sizes of and The capture 
rates can now be written as 


G ~ - C'exch(c^xx)> (27) 

C ~ Cg{(T^^) + Gexch(crxx), (28) 

by using that Geject(cr) — Gexch(o') and one hnds that D takes the form 

D = C-C = Cg{(j^^) - Cg{(j^^) - 2Gexch(axx)- (29) 

Note that even if the scattering cross sections and are equal, D will be non-zero. 
This is due to the fact that ejection of the more dominant species occurs at a larger rate. 

IV. 1. Asymmetric capture and A 

Figure 1^ shows the size of A over time for a DM mass of 5 GeV, a capture rate of DM 
on solar nuclei at c = 10^® s“^ and a capture rate of anti-DM on solar nuclei of c = 0 s“^, 
and various different self-scattering cross sections. It can be seen that, when the capture of 
anti-DM occurs primarily by DM (negative D), the asymmetry is smaller than if there would 
be no self-capture at all or the difference between and is such that D is small. On 
the other hand, the exponential growth is apparent when is larger than as to make 
D positive. Since A is dehnitely smaller or equal to A, the geometric limit of self-capture 
has dehnitely been reached once cr^x^^ > ^’"x' ^ redehnition of C and 

C would have already been necessary for the two cases with D > 0. 

It is of interest to compare the numbers captured in a Dirac DM model to those of an 
ADM model which is shown in £g. for two capture rates of anti-DM on solar nuclei. It 
can be seen that even if c = c/2, the total number of captured particles in an ADM model, 
AAsym, is at most an order of magnitude larger than the number of captured particles in a 
Dirac DM model in a large region of the plane. Note that in the case that c = 0 
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Time [Gyrs] 


FIG. 3. The evolution of the asymmetry A over time for a 5 GeV DM particle with a solar capture 
rate of c = 10^® s“^ and c = 0 and various and Blue line: = 0 and = 3-10“^'^ cm^, 

red line: = 0 and = 0, solid black line: = 3 • 10“^^ cm^ and = 0, dashed black 

line: = 4 • 10“^^ cm^ and = 0. The purple line shows the geometric self-capture limit for 

= 2 • 10“^^ cm^ and = 0 



c = c/2 



FIG. 4. The total number of captured DM and anti-DM particles of mass = 5 GeV with a 
cross section on solar nuclei ctsd = 10“'^^ cm^ in the plane normalized to the number 

captured by an ADM model with the same DM mass and cross section and twice the halo density 
{py. = 0.3 GeV/cm^). The cross section of anti-DM on solar nuclei is 0 in the left plane and crsD/2 
in the right. 
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Time [Gyrs] 

FIG. 5. The evolution of the stochastic asymmetry A over the lifetime of the Sun for various 
and using a 5 GeV DM mass and a capture rate at 10^'^ s“^. Blue line: = 10“^^ cm^ and 

^xx ~ ^xx ~ ^xx ~ solid black line: = 2 • 10“^^ cm^ and = 10“^^ cm^, 

dashed black line: = 2-10“^^ cm^ and = 0. 

and is small, Ai^Asym is slightly larger than half that of the captnred nnmbers in a Dirac 
model. This is dne to the fact that self-captnre, CN, is limited to a mnch smaller rate than 
c. The captnre of anti-DM is so small that annihilation barely occnr and the captnre is 
described approximately by an ADM model in which the backgronnd density of DM py is 
halved. 


IV.2. Symmetric capture and A 

The case of c = c implies a simple solntion to the steady state of eqs. Q and ([^. When 
the capture rates are equal, the symmetry of the equations implies that N = N. Both the 
capture of DM and anti-DM will then be given by 

N = c+{C + C)N -TN^. (30) 
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In the steady state, the capture rate is equal to the annihilation rate so that N = 0 which 
plugged into the above gives 


_C + C y/{C + C')2 + 4cr 

“ 2r ^ 2r 


(31) 


The total number of captured particles will then be 2Noo. 

Fig. I shows the evolution of A over time for various D with a DM mass of 5 GeV 
and a capture rate of 10^^ s“^. If Dt is small, the stochastic asymmetry at this point in 
time would be A = 1.7 ■ 10^^ while increasing D on the negative side will always lead to a 
stochastic asymmetry that is smaller. Considering the case of positive D, the asymmetry 
will increase exponentially with D. With the same DM mass and solar capture rate but 
turning on self-capture using = 2 ■ 10“^^ cm^ and = 10“^^ cm^, the number of 
particles at which equilibrium occurs is N^o = 2.98 ■ 10^^ while the size of A is 4.9 ■ 10^^. 
With the same self-scattering cross section and setting = 0, the stochastic asymmetry 
becomes 5 orders of magnitude larger than the expected number of particles which slightly 
decrease to N^o = 2.93 ■ 10^^. Since A is dehned as the standard deviation of A, we may in 
this case expect the asymmetry to be large. However, if the actual asymmetry was of this 


size, the geometric limit has kicked in and A is no longer given by eq. (12). Still, this is an 
indication that self-capture combined with a stochastically induced asymmetry has lead to 
a signihcant accumulation of DM in the Sun. 


V. SUMMARY AND DISCUSSION 

In this paper, we have considered the capture of DM in the Sun under the assumption 
that it is a self-interacting Dirac particle and that the galactic background density consists 
of equal amounts of DM and anti-DM. This opens up the possibility that a large difference 
in the captured amount of each type might occur so that the total number of particles in the 
Sun may continue to grow even though the annihilation cross section is that expected from 
standard thermal relic DM. The initial asymmetry between the number of captured DM 
and anti-DM particles can occur either due to different scattering cross sections for DM and 
anti-DM on solar nuclei or due to stochastic fluctuations in the capture process. Any such 
asymmetry may then be amplihed by self-capture or counter-acted by capture of anti-DM 
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by DM. The size of the asymmetry is independent of the annihilation rate and an analytical 
expression for its size was derived. 

When the capture rates of DM and anti-DM are different (c 7 ^ c), we have the case of 
asymmetric capture. If the capture rates of DM by DM is C and the capture rate of DM 
by anti-DM is (7, we define the difference in these rates as D = C — C. When D < 0, the 
capture of anti-DM is more efficient than DM self-capture which implies that the asymmetry 
A will at some point hnd an equilibrium. On the other hand, if D > 0, DM will capture 
itself more efficiently so that any initial asymmetry will grow exponentially. The asymmetry 
is found to become large enough to conclude that the geometric bound on the self-capture 
needs to be taken into account for a wide range of solar capture rates. This is due to an 
exponential dependence on the size of the DM on DM and DM on anti-DM capture rates 
when D > l/t©. This occurs for a 5 GeV DM particle with a capture rate of 10^® s“^ on 
solar nuclei, corresponding to a spin-dependent cross section of 10“^® cm^, for a D given by 
= 2 ■ 10“^'^ cm^ and = 0, where is the DM self-scattering cross section and is 
the scattering cross section for DM on anti-DM. When the capture rate of anti-DM by DM 
and the capture rate of DM by DM makes D negative, the asymmetry will always be smaller 
than in the case when there is no self-capture at all. In any case, the size of the asymmetry 
implies that the total amount of DM captured may be large. Taking the geometric limit into 
account and comparing the number of captured DM and anti-DM particles in the Sun by 
those of an ADM model, it is found that in the region of parameter space where 
a sizeable population of captured DM occurs where the numbers are within an order of 
magnitude than those of an ADM model with similar parameters but a twice as large halo 
density. Since ADM models may have an impact on solar observables as demonstrated in 
na, it is plausible that Dirac models may also have an impact on these. However, in regard 
to solving the solar composition problem and constraining Dirac DM models, numerical 
investigations such as those found in [131 El 129] needs to be performed. 

When the capture rates on solar nuclei are equal for DM and anti-DM (c = c), the 
stochastic asymmetry A, which estimates the typical magnitude of the actual asymmetry 
induced by the stochastic variation of c and c, is always extremely small in comparison to 
the total number of particles in the Sun when the self-scattering cross sections and 
are such that D is small or negative. However, in the case when D is positive, the 
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exponential dependence of D may bring the stochastic asymmetry to a size several orders of 


magnitude larger than the expected total number of trapped particles at steady-state with 


no asymmetry. However, this case is an extreme since the scattering cross section is 
taken to be right around the upper bound while = 0. Increasing to half that of 
the stochastic asymmetry is reduced by over 10 orders of magnitude to a negligible level 


compared to the expected amount for symmetric capture. The window for the asymmetry 


induced by stochastic variations for the Sun is very small and requires so it 

may not be expected that the solar asymmetry is large. However, the capture rates increase 
proportionally to the background density of DM so that larger self-capture rates may be 
expected for stars in regions where the background density is larger. Even if the likelihood 
that the Sun has a negligible asymmetry since the background density is small, stars in such 
regions may have a D that is several orders of magnitude larger which would increase A 
signihcantly thus affecting the evolution of such stars. 

In this work, we have neglected the fact that a stochastic capture rate on solar nuclei 
may imply a stochastic variation of the self-capture rates as well. One such scenario would 
occur if there are perturbations in the local background density. The investigation of this 
case is left for future work. 
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Appendix A: Capture rates of DM 

The capture of dark matter in celestial bodies is a standard calculation, hrst done by Press 
and Spergel |1], later improved and corrected by Gould [IS]. Given the velocity distribution 
of halo dark matter, /(m), in the frame of the Sun where u is the velocity very far away 
where the gravitational potential of the Sun is negligible. The capture rate is given by 



(Al) 


where wQ is the rate at which a particle with velocity w at radius r will scatter and lose 
enough energy to be captured. 
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Solar element capture 


For the SI cross section, Gould found that 

2Eo/i+ 


w^i = aini{r)- 


jj, 


u 

e 


2Bn 




(A2) 


while for the SD cross section 




wil = apnH{r){w^ - —u^)6{-^ - 


(A3) 


nj(r) is the radial distribution of element i in the Sun. The mass of and scattering cross 
section on element i is m* and <7*, respectively. The mass of the DM particle is m^, ^ 

and /i+ = The SI scattering cross section scale as 


— (^p^i 2 

H'P 


(A4) 


where Ai is the number of nucleons in element i and /Xj the reduced mass of element i and 
the DM. (7p is the proton scattering cross section and /ip the reduced mass of the DM and 
the proton. The more complex formula for SI scattering is due to the form factor which 
takes the nuclear structure of the target into account for larger energy transfers given by 




2 _ „-^EIEi 


where 


= i?. = [ 0.91 + 0.3 1 fm 


(AS) 


(A6) 


and Mj is the mass of nuclei i. For hydrogen, this form factor is set to unity. The total 
capture rate for a SI capture rate is then the sum of the capture rate by each individual 
element. For SD capture, hydrogen is the only element of importance since there is no 
enhancement of the cross sections on other elements and the fraction of elements with spin 
is completely negligible compared to hydrogen. 


Self-capture 

If DM has a non-zero scattering cross section on other DM and anti-DM particles, it 
may also be captured by colliding with other DM and anti-DM particles. A derivation of 
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self-capture is given in ref. [2l] and we review the result and add three cases for which the 
target particle is ejected. For self-capture, is broken down to 


= an{r)wPcap. 


(A7) 


Here, a is the DM self-scattering cross section, n(r) the radial distribution of already cap¬ 
tured DM and Pcap the probability that the particle is captured in a collision while not giving 
the target particle enough energy to escape the Sun. The projectile and target particles are 
gravitationally unbound when their kinetic energy is greater than m^v‘lg^/2. This means 
that, for capture of a particle without ejecting the target, the energy transfer must be 
in the interval 


u 


< 


AE 


< 


■^esc 


(AS) 


E 

The energy transfer distribution is assumed uniform on the interval which gives Q{r,w) as 


wQ = an(r)(ng^^ - u‘^)9(vesc - u) 


(A9) 


and the self-capture rate is given by 


"Rf) 


Cs = 




f{u) 


u 


an{r){y‘lg^ — v?) du dr 


(AlO) 


Ejection 

When the transferred energy in a collision involving DM and anti-DM is greater than 
^x'^Lc/2, the particle that is hit will be gravitationally unbound and escape the Sun. This 
ejection rate is calculated using the same formula as self-capture but with a different D. We 
can divide ejection into two regions, one in which u < v^sc and one in which u > v^sc- 

li u < Vesc and the incoming particle is trapped after a collision, the target particle may 
or may not be trapped. If the target is trapped, we have the case of self-capture that is 
described above. If the target particle is ejected, the transferred energy is in the range 

AE 
w'^ E 

and wQ, is given by 

wVt = an{r)u^Q{vesc — u)- (A12) 
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If M > Vesci entrapment of the incoming particle will always result in the ejection of the 
target particle. In the case that the target particle is still trapped, falls in the interval 


AE 

E 


The factor wQ is then 

wQ = an{r)vl^^Q{u - Vesc)- 

However, if AE falls in the interval 

E 


(A13) 


(A14) 


(A16) 


The target particle will be ejected and the incoming particle will still have a velocity that 
is larger than the escape velocity and thus also escape. For this case, wQ is found to be 


wf] = an{r){u^ - vl^^)Q{u - v^sc)- 


(A16) 
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